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Abstract 

An assumption of isospin symmetry permits the determination of sin(2a) from the experimental study 
of B —* tttt decays. Isospin, however, is merely an approximate symmetry; its breaking predicates a 
theoretical systematical error a 1 ^ in the extraction of a. We focus on the impact of 7r° — 77, rj mixing, as 
well as the manner in which it is amenable to empirical constraint, and determine that 0"™ can potentially 
be controlled to 0(1°). 
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I. INTRODUCTION 



Probing the mechanism of CP violation in the 5-meson system demands that the angles of the 
unitarity triangle, a,/3, and 7, be extracted a plurality of ways The study of B — > tttt decays, 
e.g., permits the determination of sin(2a), where a = (fi 2 = Arg(— V^y^/V*^/^ and we recall 
that a + (3 + 7 = tt (mod 2tt) in the Standard Model Q. The measurement of the time- dependent, 
CP-violating asymmetries in B, B — > 7r + 7r _ decays determine SVtt and C^; these measurements in 
themselves determine sin(2a e fj), where 

• /o \ 7r7r /I N 

sm(2a e ff) 



2 

7T7T 



Penguin contributions make the parameter a e s differ from a. Gronau and London have noted, 
however, that the "pollution" Aa = a C R — a can be determined and removed with additional 
B, B — ► 7T7T data under an assumption of isospin symmetry 0. Isospin is an approximate symmetry 
- the u and d quarks differ in both their charge and mass; such isospin-breaking effects make the 
determination of Aa imperfect. It is our purpose to study these effects, to the end of assessing 
the irreducible theoretical error cr™ in the determination of Aa via this method. Such is crucial to 
precision tests of the Standard Model of CP violation, realized through improved measurements at 
the current 5-meson factories and beyond |4J. Earlier work has focused on the impact of |A/| = 3/2 
electroweak penguins [^] and on the role of tt° —rj, rf mixing on Aa 0. The treatment of the latter, 
due to Gardner p], has recently been reexamined by Gronau and Zupan 0. The purpose of this 
article is to correct that work [7( and to update the analysis of Ref. 



II. ISOSPIN ANALYSIS IN B -» vrvr DECAY 

We begin by reviewing the isospin analysis in the isospin-perfect limit |3]. In this limit, the 
7r + , 7T°, and 7T _ mesons form a degenerate isospin triplet, and the B — ► 7r + 7r~ decay amplitude must 
be symmetric under the exchange of the two pions, as per the constraint of Bose symmetry. We 
can relate the two-pion states to states of definite isospin /, \(im)j), via 

I7T + 7T-) OC |(7T7r)o) + -J=|(7T7r) 2 ) 

IttV) oc |(7r7r) }-V2|(7r7r) 2 }, (2) 



where the properly symmetrized state is |/T~7r + ) sym = (l^i^) + ^1^2 ))/v / 2 = \/2|7r + 7r~) 
The particles of the tt°tt° final state are identical, so that this state need not be symmetrized - 
if symmetrization is performed nevertheless, an additional factor of 1/2 must be applied to yield 
the correct branching ratio. The B~,B° form a degenerate isospin doublet as well; evaluating the 
Clebsch- Gordon coefficients allows us to write a decomposition in terms of amplitudes Aj of definite 
isospin |8J: 

A B o_ 7T+n - = (7i + ti~\Hw\B ) = A + "^1^2 , 



2 



A B < WV o EE (7T 7T |^|5 ) EE Aq — y/2A 2 , (3) 

3 

Ab+^+tA ee (n + n°\H w \B + ) = -A 2 , 

where analogous relationships in the charge-conjugate modes (A) are implied. We recall that Aq 
and A 2 are generated by |AJ| = 1/2 and |A/| = 3/2 weak transitions, respectively. Since the 
symmetrized states |7r~7r + ) sym and |7r + 7r°) sym appear in the physical amplitudes.the B — > 7r + 7r~ 
and B — > 7r + 7r partial widths are a factor of 2 larger than suggested by Eq. © 9]. We note that 
the reduced transition rate ^{B — > 7Ti7T2) is related to the partial width — > 7ri7r2) via 



i / ^ (m^ +M 7r2 )n / _ {M m -M W2 y 
r 2 B A ^1 

Employing experimental masses throughout, we find, in specific, 



7+ _ EE 7 (B -> 7r + 7r) = 2|A B ^+,-| 2 , 

7+0 EE 7 (£? + - 7T + 7T ) = 2| A B+ ^ +7r | 2 , (5) 
7o0 EE 1 (B -> 7T 7T ) = I^B-W) |* ■ 

We note that electromagnetic radiative corrections, which may well be important, should be applied 
to yield the empirical decay widths and ultimately the reduced transition rates 10]. Irrespective of 
such corrections, if we rewrite the amplitudes of Eq. 0, which satisfy the "triangle relation," 

—=(AB-nr+ir- ~ ^B-*tv°tt q ) = ^B+^w+ttO , (6) 



V2 

in terms of the amplitudes 



>7T + 7T~ 


— A B ^ W 


A + - 


EE V2A 




= -A B 




EE V2A 



B- 



(7) 



l B+^TT+TT° , 

defined so that \Aij\ = y/jij, we find 

A+_ + V2A 00 = V2A +0 . (8) 
We assume "v/2|v4 2 | > \A \, so that i^o^o < 0. This assumption is consistent with theoretical 



assays of B — > nn decay in an operator product expansion framework [11[, as well as with the 
pattern of empirical branching ratios, given current errors [3]. Note that a similar triangle relation 
holds for the charge conjugate modes and thus 

+ V2A 00 = V2A_ , (9) 

where we note |A_o| 2 = 7(-B~ — > 7t~tt°). The form of Eqs. (JBJ) and © is identical to that of earlier 
analyses jESQ; once differing definitions are taken into account. The upshot of the isospin analysis 
is that the shift of a e ff from a induced by the penguin amplitude in B — > 7r + 7r~ , namely 

Aa ee a eS - a ee ^Arg (e 2 ^A + „A* + j , (10) 
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can be expressed in terms of empirically determined quantities. In particular, with = 

Arg(v4_| — A*i_q) and <fi = Arg(A_| — A*_ ), we have 

At>« pi , = \(i - <t>) = \ (Arg^I+^J - Aig(e*M-o,4%)) . (11) 

In the isospin-perfect limit, as we examine here, A +Q = exp(2z7)A_ 0) so that the second term 
vanishes and Aai SOS pm = Aa. Geometrically this implies that the two triangles share a common 
side; namely, t/t+o = y/l-o- Let us now turn to an analysis of isospin-breaking effects. Interestingly, 
the most significant uncertainty arises from the manner in which Aa isospin 7^ Aa, promoting the 
importance of direct theoretical assays of Aa [f^ |. 



III. ISOSPIN BREAKING IN B -» tttt DECAY 



The charge and mass of the up and down quarks do differ, so that the predictions of the isospin 
analysis we have discussed cannot strictly hold. There are two different effects to consider. Firstly, 
penguin contributions of |AJ| = 3/2 character can occur, mediated either by electroweak penguin 
effects, or by isospin-breaking in the strong-penguin matrix elements 0, Q, 0]. Secondly, the 
triangle relationships of Eqs. (jH)) and © need no longer hold For example, the physical, 
neutral-pion state contains isoscalar components due to mixing with the 77 and 7/, engendering an 
"J = 1" amplitude in B — > tttt decay p. The 77 and rf admixtures in the tt° are generated by 
the strong interaction in Oim^ — m u ). Alternatively, we can regard this interaction as an / = 1 
"spurion" , encoding isospin-violating effects so that the matrix elements with the spurion are SU(2) f 
invariant [16(. In this latter picture, the "extra" amplitude engendered by n° — rj,T]' mixing can 
be recast as a |A/| = 5/2 amplitude, generated by 0{md — m u ) or 0(a) effects in concert with a 
I AI I = 3/2 weak transition. A |A/| =5/2 transition can also be realized by isospin-breaking effects 
in concert with a |A/| = 1/2 weak transition |9j, though, as per the spurion picture, such is not 
engendered bv 7r° — mixing in leading order in isospin breaking. Writing A|a/|,/, we replace 



Eq. © with $ III 



A B 0^>TtH° 
Ab+^tt+ttO 



Al/2,0 + —=(^3/2,2 - ^5/2,2) , 
Al/ 2 ,o - ^(^3/2,2 - A 



15/2,2 J 



3 . 3 . 

2 A 3/2,2 + y 2^5/2,2 



(12) 



We note that this parametrization suffices to capture isospin breaking in B — > tttt decay, as three 
theoretical amplitudes describe the three empirical ones. Isospin breaking impacts the determination 
of Aa in two distinct ways. For example, it can break the triangle relation, Eq. ©. If the triangle 
relation is broken in an ill-determined way, the ability to assess the angles (f) and (f) is compromised. 
If the triangle relation is not broken, however, then the application of the isospin decomposition 
given in Eq. (jSJ) permits the determination of and regardless of whether additional isospin- 
breaking effects are present. Isospin breaking, however, can also make Aai SOSp in differ from Aa; 
specifically, Arg(e 2 * 7 A_ A^ ) 7^ 0, recalling Eq. (fTTj). If the impact of both effects can be estimated, 
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if not controlled via empirical constraints, we can assess the irreducible theoretical error in the 
determination of Act. Interpreting these two effects in terms of the parametrization of Eq. (|12jl. a 
non-zero value of the A 5 / 2 ,2 amplitude signals the breaking of the triangle relation, whereas penguin 
contributions to Ab+^+^o, to A 3 / 2) 2, make Arg(e 2l7 A_oAl a ) 7^ even if A 5 / 2 ,2 = 0. Electroweak 
penguin contributions are an example of the latter effect |l8j]. Since current experimental data is 
consistent with | ^.3/2,2 1 ^ I ^1/2,0 1 m B — > nn decay, we expect that tt° — r},rf mixing will play the 
most important role in the realization of a A 5 / 2i 2 amplitude. In the limit that the A 5 / 2j 2 amplitude 
is generated exclusively in this manner, Arg(e 2r) ' A-qA* +q ) 7^ can only be realized through penguin 
contributions of |AJ| =3/2 character. The phenomenon of 7r° — 77,77' mixing can generate both 
effects; let us consider it explicitly. 



A. 7r — 77, 7/ Mixing 

In what follows we examine the role of tt° — 77, 77' mixing on the extraction of a from B — > ttti 
decays. We distinguish the amplitude for decay to physical pion final states, which suffer 7r° — 77, 77' 
mixing, e.g., ^^0^0, from the amplitude in the isospin-perfect limit, Ag^^, where 3 denotes the 
isospin-triplet state with I 3 = 0. Noting earlier work on 7r°— 77, 7/ mixing in K — > inr decay 11130, 
we have 

A B+ 

— >7r+7r° — Ab+^tt+^s 4~ eAb+^jt+t) ~l~ £ Ab+^-k+t] 1 5 

Ab^tt tt = AB^fofo + 2eAB^(j, 3V + 2e' Ab^^ti' , (13) 

where we assert e,e' ~ 0((rrid — tti u ) / A^) or 0(a) and neglect all higher-order terms in isospin- 
breaking parameters. To gain insight on the nature of Ah a d ; we note that the analysis of the 
pseudoscalar meson octet in current algebra [ijj, or in lowest order chiral perturbation theory |22| |. 
determine the ir° — r]s mixing angle eg to be 

/ m d - m u \ 

e 8 = — ( , 14 

4 \ m s — m J 

with 777 = {m u + m d )/2, so that we expect A had ~ 0{m s ). The impact of isospin breaking is 
controlled by the magnitude of SU(3) / breaking. The breaking of SU(3) / symmetry also engenders 
the mixing of the pseudoscalar octet and singlet states, rj 8 and 770, to realize the observed 77 and 77' 
states. Such considerations demand that we evaluate A B ^ nv (i) in the presence of SU(3)/ breaking 
effects. We postpone specific estimates of e and e' to the discussion of our numerical results. We 
may use these relationships to rewrite the triangle relation, Eq. ©, which now appears as 

-^(Ab^tt+tt- — Ab->4, 3 4, 3 ) = A B +^TT+(j,3 , (15) 

in terms of amplitudes employing physical 7r° states. That is, 

-^(Ab^tt+tt- — A B ^o n o) = A B +^+ n o — V2eA B ^ 3V — V^e ' A B ^ sr f 

—eAB+-* 7T + v — e As+^TT+rf , (16) 
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where replacing Ab^^-co with A B _^o„w generates corrections of higher order in e,e', which are 
negligible in the order to which we work. Note that 77^ connotes either 77 or 77' throughout. We 
observe that the triangle relation is broken in the presence of isospin- breaking effects [6J. 

We turn to theory to assess the impact of the amplitudes containing 77, 77' on Eq. (fTKj) . The 
QCD factorization approach (23| to hadronic 5-meson decay analyzes the decay amplitudes in a 
systematic expansion in inverse powers of the heavy-quark mass m& and the strong coupling constant 
a s (fi), where \i ~ 0{m h ). Crucial to the treatment of the decay amplitudes in this case is that of 
the physical 77 and 77' states themselves, as the 77 and 77' mix. These states are not simple flavor-octet 
and flavor-singlet states, as SU(3)/ symmetry would suggest, but rather each physical state is a 
mixture of these components. The presence of the flavor-singlet component in decays to 77 W final 
states admits novel decay mechanisms, in part mediated by the axial anomaly, not present in other 
channels j^ij. We emphasize, as recognized in Ref. j^ij, that these flavor-singlet contributions are 
not captured by the analysis of non-77 ^ decay channels with an assumption of SU(3)/ symmetry. 
To implement 77 — 77' mixing, we employ the Feldmann-Kroll-Stech scheme J26|, also adopted in 
Refs. J24], EEj|, in which a single mixing angle characterizes the decomposition of \rj^) into the flavor 
states \r] q ) = i\uu) + \dd))/\/2 and |?7 S ) = \ss). Beneke and Neubert thus determine ^T* 



V2A™ 
+V2A< 



, 94 BN 



p=u,c 



mis 



3 3 



+A w 

Vq ' 



3 3 3 3 

Spuiat + a 2 + Pi + fa) + 7> a l,EW + o a l,E\v + o@Iew + ~PIew 



(17) 



where A r 



V p bV* d , and we refer to Ref. [25] for all details. In this expression, the meson masses, 



decay constants, form factors, and light-cone distribution functions are all evaluated in the isospin- 
symmetric limit. Note, however, that the physical quark charges have been employed, so that, in 
particular, e u 7^ e^, in the evaluation of the electroweak penguin contributions. We note that the role 
of a possible cc component in the 77 ^ mesons has been included in the computation of Ar-^- m 
and Ago_ H 4 3 n(o, although such effects are likely most significant in b — > sqq transitions [271] . which 
we do not treat here. The amplitudes are computed in next-to-leading order in ot s and at leading 
power in AQCD/ m 6- The terms representing weak annihilation contributions are denoted by and 
are included, although they are formally suppressed by a power of m b . The computation of these 
contributions suffer endpoint divergences in QCD factorization, so that their estimate is uncertain. 
The large direct CP asymmetry found in the penguin-dominated mode B — > K + ti~ sug gest s 
that annihilation contributions may well play a larger phenomenological role than anticipated [25]. 
Nevertheless, if we do neglect the annihilation contributions, as they are power-suppressed and 
largely possess, in this case, the same weak phase as the dominant contributions, we have 



V2Af__ 



>7T 77 



(') 



+ 2A 



BN 



E 



p=u,c 



A, 



A co 

Vq 1 



-a 



p 

3,EW 



+ —a 



p 

A,EW 



'11 



where the a, implicitly depend on the order of the arguments of the Aij prefactor, so that a\ 



(p) 
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a^(r/g V). By comparison, we note that 



V2AT 



A m 



(19) 



jpuy^l l ^A) i 2 3,-EW ' 2 

p=u,c % 

where we emphasize that «; p) =«! p) (7r7r) in this case. The electroweak penguin contributions which 
appear in this expression are explicitly of |AJ| =3/2 character. That is, were the quark-charge 
dependence made manifest, we would see that these contributions are proportional to e u — e^, so 
that, in analogy to our discussion of tt° — T),rf mixing, the electroweak contribution contains an 
effective isovector interaction acting in concert with a |A/| = 1/2 transition. Thus we can write 



.4 



BN 



A BN y 



X (/) 

Vq 



A (,) 
Vq * 

A 



E(7T7r) 

3 3 
5 pu (ai + a 2 ) + -al EW + -a p 4EW 



(20) 
(21) 
(22) 

satisfy both Eqs. © and (JSJ), though the Aj thus determined would be v^2 larger. Nevertheless, no 
physics can depend on this normalization choice, so that we must have 



where 



and 



E(MiM 2 ) = ^2 \ 



p=u,c 



with af ] = af ; (MiM 2 ). We note 



» 



A^j?_^ j k /y/2, so that the amplitudes AW 



I ABN I 



1-4 



B--+7T-TjW I 



1-4 



_B — >7T </>3 | 



I yjBJV I 

I ABN f 

I B°-k/> 3 (/>3 I 



1-4 



B°^</> 3 ??(') 



1-4 



B°^0 3 03 I 



as well as 

Ab--+ti-»j 

Returning to Eq. (JTfij) . we find 

1 



B -K/> 3 r)<'> 



v/2 



(4b— tt+tt" - ^B-»7r 7r°) 



BH 



where £, which need not be real, is given by 

Z = e[X Vq + ...]+e' [X v ,+...] 



(23) 



(24) 



(25) 



(26) 



— (jp — Qp Qp 

We note that X (/) < — ► X w, £(MiM 2 ) « — > £(MiM 2 ), and A„ < — > A* under CP transformation, 

Vq Vq " 

whereas each ellipsis denotes neglected annihilation corrections. It is worth emphasizing that this 
result differs from its analogue in Ref. Q in an important way. That is, we have not assumed 
SU(3)/ symmetry in the construction of £, whereas Ref. 0] neglects all SU(3)/-breaking effects save 
for 7] — 7]' mixing. In Ref. Q, £ is replaced by the parameter eo, namely 



(27) 
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Let us examine the ingredients of Eqs. (|26j) and (|21|). The ratios of S(MiM2) differ from unity if 
SU(3)/ is broken in the light-cone distribution functions, through, specifically, the hard-spectator 
contributions to ai(MiM 2 ) The latter are real if the contribution of the twist-3 distribution am- 
plitudes, which generate divergent, albeit formally power-suppressed, contributions, are neglected. 
This implies that the ratios of S(MiM2) in Eq. (}2*T]) are real, if all subleading corrections and 
electroweak penguin effects are neglected. The latter do make the ratios of S(M!M 2 ) complex in 
leading order in l/m&; however, electroweak penguin effects enter £ in 0(a{m d — m u )/A had ), so that 
their inclusion is actually of higher order in isospin breaking. The remaining factors in Eq. (}2*T|) are 
given by 



A (,) 

A 



iB- 



^ ( "(0) 



F B- 



-(0) 



(28) 



where F S 



denotes a form factor for the decay to a pseudoscalar meson M in the convention of 



Bauer, Stech, and Wirbel 28] . Following Beneke and Neubert, we parametrize [24 



(29) 



noting that F\/F% ~ 0(1) in the heavy-quark limit. The first term is related via SU(3)/ breaking 
to F(f^(0), where one expects F\ ~ F^ n {fi) in the FKS scheme jiij]. The second term, however, 
is driven exclusively by the flavor-singlet contribution and cannot be related to -F ( ^ _>7r (0). It is 
ill-known, though likely of greater impact on the B — > rj' form factor 



1. Breaking the triangle relation 



A non-zero value of the parameter £ in Eq. (j2*3j) signals the breaking of the triangle relation, 
Eq. ([B|L and the appearance of an amplitude of |AJ| = 5/2 in character. In the QCD factorization 
approach, £ is given by Eq. ([2l>]). If £ can be determined with surety, and is real, an "isospin 
analysis" based on Eq. (j23j) can determine <fi and without theoretical error from this effect. 
Determining £ requires the isospin-breaking parameters e and e', which characterize tt° — 7], rj' 
mixing, as well as X w, as per Eq. (12111. We begin by determining X w in the QCD factorization 

Vq ~ * Vq 

approach. The parameter X^o is controlled by F^ 11 / Fq^ 77 and S(?7q7r)/S(7r7r) exclusively, if 
power-suppressed contributions are indeed negligible. The former drives the numerical value of £. 



Using the parameters of Ref . [24J , we find 



-(0) 



F B- 



(0) 

F^-(0) 



0.83 ±0.02 -> 0.89 



0.68 ±0.02 -> 1.1 



0.82 



0.58 



(30) 



where the reported errors are determined from the errors in the inputs alone, assuming they are 
uncorrelated. The first number reported for each ratio employs F 2 = 0, whereas the second number 
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employs, rather arbitrarily, F2 = 0.1 as per Ref. |24|. For reference, we have also included, in 
brackets, the ratios assumed in the SU(3Vf approach of Ref. Q. The form factor ratio for the 7/ 
can differ substantially from that of Ref. jjj , and varying the 77 — 7/ mixing angle does not capture 
the excursion found. As for the remaining factor, we estimate, neglecting power corrections and 
electroweak penguin contributions, 

S(r/ g 7r) ^ ai(^V) + Q^qV) 

S(7T7r) ai(7T7r) + a2(7T7r) 

* 1 + i^r 1 + ^)(«? - «J) (si) 



1-5-10 



-3 



where the deviation from unity is determined by SU(3)/ breaking in the light-cone distribution 
functions, as parametrized by alf, which appear in the hard spectator terms. We note, as in the 
case of annihilation contributions, that endpoint divergences can appear in the power corrections. 
We employ the parameters given in Ref. |25( and observe that this source of SU(3)/ breaking 
appears to be negligible. This observation is consistent with other recent data. For example, 
SU(3) f breaking in the form factors and decay constants suffices to ex plai n the large difference in 
the observed branching ratios for B s — > K + K~ and Bd — > tt + tt~ decays We thus determine 

X Vq = 0.83 ±0.02 -> 0.89, 

X v , q = 0.68 ±0.02 -> 1.1. (32) 

It is worth noting that the form of Eq. ((23$ is quite general; it does not rely on our adopted 
framework for 77 — 77' mixing. If, instead, a general, two-angle mixing formalism |3lL H3| in the 
octet-singlet basis were employed to describe 77 — 77' mixing, an equation of form Eq. (J21j) would 
nevertheless emerge Jiij]. We would also find compatible numerical results 1 . Indeed, we can use 
the empirical decay amplitudes to define and determine an effective parameter X c f r] via 

Vq 

AB--*n-T)W + V2A s ^o v (o = A B --^- ir oX^ ) , (34) 
where if power corrections, as well as isospin-breaking effects, are negligible, X e f n is X m as de- 

Vq Vq 

fined in Eq. f!21|) . Empirical branching ratios for B~ -> vr-r/W, B -> tt^W, and Br — > IT 7T 



1 In a two-angle mixing formalism, assuming as per Ref. j3| 1 that the ratios of the B — > jjW and B — > n form factors 
are determined by the ratios of the related decay constants, we have 

f B ""(0) = V2/ f s cos9 8 /osingp ^ 
F B ^'(0) _ V2f f 8 sm0 8 /qCos^ 



If /g = /o = /tt and 0s = 0q = 9 and we assume the ideal mixing angle = sin x ( — 1/3) we recover 
F o~* n (°)/ F Q^(°) = VV3 and F ( f^ r ''(0)/ J F s ^ 7r (0) = JT/3 as per Ref. 7]. If we employ the parameters 
in either Eq.(3.7) or Eq.(3.8) of Ref. |22| we find results comparable to what we have reported in the F2 = case. 
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TABLE I: CP-averaged branching ratios for selected B — > PP modes from the compilation of Ref. [l2j], 
reported as 10 6 Br(B — ► PP). We display the experimental data, both preliminary and published, available 
since the compilation of Ref. [3^ ] and included in the averages of Ref. . 



Mode 
PP 


PDG [36] 


BABAR 


Belle 


CLEO 


CDF [37] 


HFAG [12] 


TT+TT 


5-6±?;? 


5.8 ±0.6 ±0.4 [38] 


5.0 ± 1.2 ±0.5 [39] 


4.6±^t° ; 6 r [40] 




5.5 ±0.6 




< 5.7 


5.1 ±0.6 ±0.3 [41] 


4.8 ±0.7 ±0.3 [42] 


1.2±f;| [43] 




4.9 ±0.5 




< 7 


4.0 ±0.8 ±0.4 [41] 


< 7 [44] 


i.o±5:g [43] 




4.0 ±0.9 


7T + 7T~ 


4.8 ±0.5 


4.7±0.6±0.2 [45] 


4.4 ±0.6 ±0.3 [39] 


4.5±1;1±«;| [40] 


4.4 ± 1.3 


4.5 ±0.4 


7T°7T° 


1.9 ±0.5 


1.17 ±0.32 ±0.10 [38] 




< 4.4 [40] 




1.45 ± 0.29 


W ° 


< 2.9 


< 2.5 [47] 


< 2.5 [42] 


< 2.9 [43] 




< 2.5 


77'tt 


< 5.7 


< 3.7 [4JJ 




< 5.7 [43] 




< 3.7 



decays can thus eventually determine l-X" 6 ^!, and the angles AYg(A B -_ i7T - rj (,)A* B __ t7T _ 7T0 X e ( ,f) and 

Vq ' Vq 

Arg(A^__ >7r o^(/) A* B _^ 7T _ w0 X e (l f), up to discrete ambiguities. Data on the charge conjugate modes 

Vq 

would determine the charge-conjugates of these quantities in a similar manner. Our theoretical 
analysis suggests that X e f is real to a good approximation 2 , so that the deduced empirical angles 

can be interpreted as Arg(y4 B -_^ 7r -^(,)y4^__^ 7r _ 7r0 ) and ATg(A B ^o v (>) A B _^_^ ) 0- Nevertheless, 
verifying that IX ^! = IXtjJ, e.g., would serve as a consistency check. We thus expect that this 

Vq Vq 

analysis would not only constrain the ill-known B — > rf 1 ' form factors, but also help determine the 
extent to which Ao;i Soa pi n 7^ Act, as we shall explain. 

Before turning to this issue, let us conclude by determining the expected value of the |A/| = 5/2 
parameter £ and the manner in which its uncertainty impacts Aai SOS pin- To compute £, we use the 
recent results of Kroll for the 7r° — 77, 7/ mixing angles [34j : 

e = 0.017 ±0.003 ; e' = 0.004 ± 0.001 , (35) 

to yield 

f = 0.017 ±0.003 -> 0.020 [0.016], (36) 

where the error in £ is determined from those of the inputs alone, assuming their errors are uncor- 
rected. We have incorporated the ir° — 77, 7/ mixing angles directly as determined in low-energy 
experiments; we note that the scale dependence of the light-cone distribution functions, of which 
this is part, does not enter at next-to-leading order accuracy in a s [23]. The results employ F 2 as 
in Eq. fl3TI|) and report, in brackets, the value found in the SU(3)/ analysis of Ref. 0] as well. Note 
that the greater uncertainty in X v > noted previously has little bearing on the final error in £, as 
the 7T° — 77' mixing angle e' is relatively small. Given an estimate of £ and its error, we can also 



2 The SU(3)/-breaking electroweak penguin contribution to £(77,7 )/T,(tttt) generates the only complex contribution 
in leading power in 1/mi,. Note that of the neglected annihilation terms, only the electroweak penguin annhilation 
contributions can be complex. 



10 



proceed to determine the uncertainty in consequent to it. To do this, we note that cos0 can be 
determined from the empirical decay amplitudes, determined from the empirical branching ratios 
via Eqs. (|4I5|) . and the relationship given in Eq. (J2H1) 3 - We employ Eq. (JHJ) for the neutral modes, 
but define in this case 



:i-oi^ 



to yield 



where 



cos< 



cos< 



= cos q 


>o + £ 


( cos O V2 '^^'V 


- o(e) , 


1 


\\Ab- 


.7T + 7T- 2 — A B ^, n n 2 + 2 A B + 


12 

-^7r+7T | 


2^2 









(37) 



(3* 



(39) 



The error in £ generates an error in the determination of 0; namely, a§ = a^\d<f)/d^\, where we 
note that the error in follows from replacing the amplitudes by their CP conjugates. Assuming 
a 100% error in our estimate of <7g, as F2 is ill-known and the ir° — r),rf mixing angles can have a 
small electromagnetic component, estimated to be some 6% of £s |2pj ]. we employ 0^ = 0.006 and a 
recent empirical compilation of CP-averaged branching ratios [12| , reported in Table HI to estimate 
0$ = 0.4° and thus an error in Aai SOSp i n of 0.4°, as we add the errors linearly. In constrast, the shift 
in due to the contribution is 1.2°. We note, in particular, that our error estimate can be 
made more robust, if not reduced, through the measurement of B — » nrf 1 ' decays. 




T 



FIG. 1: Schematic illustration of the triangle relation in B — > tttt decay in the presence of isospin breaking, 
Eq. (|25jl. Note that A = A B _, n + n - / ^/2, B = A B _, n o 7T o /\^2, and C = (1 — ^)A B + _^^+^o. Moreover, £ is real, 
with A B +^ 7T + 7T o = T + P, where T is the tree-level contribution to B + — > 7r + 7r° decay in the isospin-perfect 
limit. 



2. Breaking Aoi sosp i n = Aa 

Thus far we have determined the error in Aai SOSp i n incurred through the uncertainty in the 
parameter £. The error in Aa, however, is determined by that in 

Aa = A« isospin + iArg(e 2 ^_ A; ) . (40) 



The sign of cf> is undetermined. 
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Defining ( = Arg(A B +_> n + n oT*) and ( = Aig(A B -^ n - n oT*), where T is the tree-level contribution 
to B + — > 7T + 7r° decay in the isospin-perfect limit, so that |T| = |T|, we have 

Arg(e 2 ^A_ A; ) = C - C • (41) 

The interplay of these relationships is illustrated in Fig. ^ As in the isospin-perfect case, the angle 
(f) has a discrete ambiguity, as does the angle (: their overall sign is not determined. This is realized 
as an ambiguity in the orientation of AABC, that is, whether it points up or down. A similar 
ambiguity also exists for the charge-conjugate amplitudes, which yield and (, to yield a four-fold 
ambiguity in Aa. The angles ( and ( are non-zero if penguin contributions of |A/| = 3/2 character 
occur. We wish to estimate the extent to which Aa ^ Aa isospin , as well as its uncertainty, though 
we shall begin by considering the contribution from 7T° — rj, rj' mixing exclusively. 

In the presence of 7r° — rj, 7/ mixing, we have 

Ab+^>tt+tt = As+^+fo + eA B +^ 7T + rj + e A B +-rK+rf , (42) 

with A B +^ n + ( f >3 = T + P ew , where we emphasize that the amplitude computed with the 1=1, 
Is = state 03 can contain an |A/| = 3/2 electroweak penguin contribution, P ew , in addition to a 
tree-level contribution T. The angle ( can be written as 

C = Arg(A fl+ _ w+ff oA^ + ^ r+03 ) + Arg^B+^+^T*) , (43) 

where the first term is of 0(e, e') and the second is rendered non-zero by P cw . Here we focus on the 
first term, namely, 



Cr h rf = Arg(A B+ ^ 7r+7T oA B+ ^ 7T+(f>3 ) = e sin 6 V g+ ^ 7r+T? + £ ' s i n 6^ 




-A— t — ^ + e sin e v > 

j(B+ -> 7T + 7T 17 




1-4 



B+^7r+r;' 



U 



(44) 



where we work in 0(e,e') throughout, with 

e v (,) =ATg(A B+ ^ 7r+v (,)A B+ ^ +7r0 ). (45) 

Defining analogous variables for the CP-conjugate amplitudes, we determine that the contribution 
to Aa — Act isospill from tt° — r/, rj' mixing is 



(Qr/ - Cvrf) — o I Sil1 9*1' 



2 



7(5"- 


-»• 7T 77) 


7 (5-- 


-> 7r _ 7T ) 




— > 7T — 77') 


/ 7(5" ^ 


— > 7T~7r°) 



7 (5+- 


> vr+r/) \ 


7 (S+^ 


• 7T + 7T ) ] 


/ 7(5+ • 


— > 7T + 7/) 


/ 7(5+ - 


-> 7T+7T ) 



Letting sin = — sin 6L(o = 1 and employing the empirical, CP-averaged branching ratios of the 
compilation of Ref . , we estimate 

\{C^-Cr,, v ')<lA ±{).2\ (47) 
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where the error follows from the errors in the branching ratios and 7r° — rj, rj' mixing angles alone, 
assuming such are uncorrelated. Note that we have added the errors in C, v>r) i and £ VtV ' linearly. 
This numerical result can be compared to the bound of 1.6° at 90% confidence level (CL) reported 
in Ref. [jj; our bound is slightly smaller as it employs measured branching ratios, rather than 
experimental bounds. Its provenance is also different, as the bound in our case does not depend on 
an assertion of SU(3)/ symmetry. Most importantly, it is improvable, as it is driven by empirical 
errors. In addition, the angles 9 „(o and 6L(,) are also subject to empirical constraint, so that a direct 
assessment of (( v ,n' ~ (W)/^ should eventually prove possible. Ultimately, it is the uncertainty in 
(Cvtf ~ C?7,?7')/2 which matters, not its gross deviation from zero. 



B. Other Isospin-Breaking Effects 

Thus far we have considered the disparate roles of tt° — rj, rj' mixing: this isospin-breaking effect 
can not only engender a |AJ| = 5/2 amplitude, breaking the triangle relation of Eq. (JUJ), but also 
generate a |AJ| =3/2 penguin amplitude, forcing Aa isoS p in — A a ^ 0. Yet isospin breaking is not 
limited to 7r° — rj, rj' mixing, and we can ask what other effects might enter, as well as how we might 
discern their presence from the experimental data. 

As we have mentioned, |AJ| = 3/2 electroweak penguin contributions can contribute to 
Aajsospin — Aa 7^ 0, through the second term of Eq. (pf3~)). If one neglects the electroweak pen- 
guin operators associated with small Wilson coefficients, namely c-j and c$ 0], then the impact of 
these contributions on Aa can be assessed, without theoretical ambiguity, up to isospin-violating 
corrections 0|, to yield 

(Aa - Aa isospin ) ewp = 1.5° ± 0.3° , (48) 

where the error arises from that in the empirical inputs. Isospin-breaking in the matrix elements 
of the strong penguin operators can also engender a |A/| =3/2 contribution , not captured 

by 7T° — r/,//' mixing. For example, corrections of 0(a) can distinguish from A^ 37T ±, or, 

specifically, j? r p ± ^ 7r± ( - g) f , from Fq ^ s (0)f n ±. In addition, m d 7^ m u effects beyond n° — 77,77' 
mixing can also occur |l4l . [l5[ . though in K — > ixtx decay, e.g., such terms do not appear in the weak 
chiral Lagrangian in 0(p 2 ) [35|. The contributions from electroweak penguin operators should yield 
the largest effect 0. In particular, we note, using the notation of Ref. [iil ]. that or | C4 1 / 1 C 9 1 ~ 4%. 

Contributions to the |A/| = 5/2 amplitude, not mediated by 7r° — 77,77' mixing, can also occur. 
For example, 0(a) effects in the evaluation of F^~™ and can yield an effective |A/| = 5/2 
amplitude from either |A/| = 3/2 or |A/| = 1/2 weak transition operators. Contributions built on 
the former can be absorbed by modifying £ and £ and enlarging their errors. Contributions built on 
the latter are more problematic, as they will make £ and £ complex, as well as £ ^ £. Consequently, 
the angle determined from the analysis of Eq. (|23j). i.e., 0' = Arg(A + _(l — £,*)A* B+ ^ w+w0 ), is not 
0. A similar conclusion emerges from the study of the charge conjugate amplitudes, where we note 
(f)' = Arg(A + _(l — £*)^4jg+_ >7r + 7r o) is not (p. Generally we can rewrite Eq. (|4T)|) as 

Aa = \{$- 0') + ~(C - + ~ [(0 - 0) - i$ - 00] , (49) 
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where the last term vanishes if £ and £ are real. Note that the geometric interpretation of 4> 
and ( illustrated in Fig. ^ as well as of <fi and (, make the signed contributions of (<f) — <f>)/2 
and (C — C)/2, and, by inference, (</>' — <$>')/ 2 and (( — ()/2, add constructively. The sign of the 
term \(<b — <j)) — (<// — <p')]/2, however, is unclear. Nevertheless, in contradistinction to K — > Tin 
decays (si EM 0| , we do expect the role of the |AJ| = 1/2 weak transition in generating an effective 
| A/ 1 = 5/2 amplitude in B — > titt decays to be a relatively small effect. That is, on general grounds, 
the pattern of empirical branching ratios shows that the |A/| = 1/2 amplitude is not dominant, 
indeed that | ^3/2,2 1 > |Al/2,2|, and a/e ~ 0.4. It is worth noting, though, that £ and £ can be 
complex from the inclusion of ir° — 77, 7/ effects alone. However, such effects arise, in leading power, 
from SU(3) / breaking in the light-cone distributions functions of the 77 ^ and tc which appear in the 
electroweak penguin contributions and, in the power corrections, through the electroweak penguin 
annihilation contributions. We find the leading-power effect to be negligibly small, though testing, 
as per Eq. (13*41) . whether \X e % I = \X e f f) I is bourne out by experiment should reveal the presence of 

unexpectedly large complex contributions. 



IV. SUMMARY 



The study of B — > txtx decays under an assumption of isospin symmetry permits the extraction 
of the angle a, modulo discrete ambiguities. This is realized through the determination of the 
penguin pollution Aa, which is also discretely ambiguous, yielding a = a e R — Aa from the directly 
measured quantity sin(2a eff ). Isospin symmetry is broken in nature, as the up and down quarks 
differ in both their mass and charge, and it is important to assess the error thus incurred on Aa. 
We have studied isospin-breaking effects in B — > mi decays, placing particular emphasis on the role 
of 7T° — 77, 7/ mixing, as it yields the most significant effects. In particular, tt° — rj, rj' mixing can not 
only engender a |A/| = 5/2 amplitude, breaking the triangle relation of Eq. (jHJ), but also generate 
a I AI I = 3/2 penguin amplitude, forcing Aa isospin — Aa 7^ 0. 

We recognize that, in nature, all flavor symmetries are approximate, and we have computed the 
shift in Aa to leading order in isospin breaking, with an assessment of the error in this shift. To 
realize this, we have worked within the QCD factorization framework, though our results do not 
depend on the details of such an analysis. Rather, the essential point is the utility of a combined 
heavy-quark, l/m?,, and ot s expansion of the theoretical decay amplitudes. We use it to sort through 
the various effects, to determine that the empirical B — ► tttt amplitudes satisfy a modified triangle 
relation, Eq. (J25|) . with a isospin-breaking parameter £ which is real, to good approximation. More- 
over, under the assumption that £ is real and determined exclusively by 7r° — 77, 7/ mixing, its value 
can be determined from experiment, once information on the 7r° — rj, rf mixing angles is employed. 
Indeed, the essential improvements over the analysis of Ref. are these: that a relationship of form 
Eq. (|23|) exists with a real parameter £ and that empirical information on B — ► irr]^ decays exists 
and can be employed to constrain the impact of isospin-breaking effects. This is important, as the 
B — > 77W form factors contain contributions which are not constrained by SU(3)/ symmetry jjl . 
The empirical data on B — > nr/^ decays is incomplete, though it can be expected to improve. 
Nevertheless, enough information currently exists to realize a crucial shift in our perception of 
isospin-breaking effects. What matters is not the shift in Aa per se, but rather the surety with 
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which we can assess that shift. We assess that the change in Aa due to isospin-breaking effects, 
namely 5(Aa) = Aa — Aao, where Aao represents the penguin pollution in the isospin-perfect 
limit, is 

6{Aa) = 1.2° [£] + 1.5° [P ew ] + 1.1° [P^_^) + . . . « 4° , (50) 

where we have resolved the discrete ambiguity in Aa by assuming that > and < 0. We 
note, currently, that a = (lOllg 6 ) no corrections from isospin-breaking effects have been 



included. The contribution labelled is the shift in Aa due to the presence of a A 5 / 2 ^ amplitude, 
realized through 7r°— rj, rj' mixing only. The contributions labelled u P ew " and u P n o_ r) ^ ri " represent the 
shift in Aa due to penguin contributions of effective |AJ| = 3/2 character. We emphasize that the 
latter number is a bound, rather than an explicit estimate. The ellipsis includes neglected isospin- 
breaking contributions, such as A 5 / 2: 2 contributions generated by 0(a) effects on the |A/| = 1/2 
weak transition, which should be rather smaller than the estimate labelled by £. However, the errors 
in these estimates are smaller and are insensitive to the manner in which the discrete ambiguity in 
Aa is resolved: 

of = 0.4° [£] + 0.3° [P ew ] + 0.2° [P^rf] + 1.1° [bound] + . . . « 2° , (51) 

and it is improvable. Note, in particular, the error associated with the v4 5 / 2 ,2 amplitude comes 
from doubling the error in the theoretical computation of £; here we employ the theoretical range 
in the F^ v() form factors recommended by Ref. j^]. This error can be tested, if not mitigated, 
through the use of anticipated empirical data. Note, too, that we have included the bound from 
penguin contributions to — > tt^tt decays from tt° — r], r( mixing in our theory error. This can 
be mitigated with improved empirical data. Ultimately, we believe that a theoretical systematic 
error cr!^ of 0{1°) is attainable. 
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